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METHOD AND APPARATUS FOR CHARACTERISING MULTIPHASE FLUID 

MIXTURES 

TECHNICAL FIELD 

This Invention relates to a method and apparatus 
5 for characterising multiphase fluid mixtures (e . g . 

slurries, emulsions, suspension of bubbles and fine solids 
in liquid, and bubble froth phase) based on the analysis 
of electrical impedance spectrum using predictive 
mathematical algorithms, such as artificial neuron 

10 network. 

BACKGROUND OF THE INVENTION 

In many industrial processes involving multiphase 
fluid mixtures where the components and mixtures may be 
stationary, moving in bathes or flowing continuously, 

15 there are needs for accurate and inexpensive phase 

concentration monitoring methods and means. It is also 
often desirable that these methods and means have the 
capability of working on-line v/ith the processes. 

A number of methods have been used in the past to 

20 monitor the phase concentration of multiphase fluid 

mixtures. Generally these methods seek to find a specific 
property which is significantly different for the phases. 
The value of this property for the mixture will then 
depend on the phase concentration. By measuring this 

25 property one would be able to find the phase 

concentration. Examples of the specific property are 
electrical properties (i.e. conductivity or capacitance), 
density, viscosity, absorption of light or absorption of 
radioactive radiation • 

3 0 Precise and safe measurement of conductivity or 

capacitance requires relatively simple instrxamentation. 
Thus, methods based on conductivity or capacitance have 
been widely used in practice for measuring phase 
concentrations not only in solids -liquid systems but also 

35 in gas-liquid, liquid-liquid and three-phase systems. 

Examples of conductivity or capacitance based 
devices are disclosed in U.S. Pat. No. 4,266,425 to 
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Allport et al., U.S. Pat. No. 3,523,245 to Love et al . 

The prior art systems described above, however, 
have a few of major drawbacks. Electrical conductivity 
based methods are very sensitive to the variations in the 
5 electrical conductivity of the liquid phase o,f the 

multiphase fluid mixture. For example, the electrical 
conductivity of an aqueous slurry may increase by more 
than 50 times with the addition of 2.5% by weight of salt 
(NaCl) to the aqueous phase. When the conductivity of the 

10 liquid phase changes siibstantially with time, the 

conductivities of both the slurry mixture and the liquid 
phase are required in order to calculate the solids 
concentration. But the on-line measurement of the 
conductivity of the liquid phase in a slurry mixture is 

15 generally difficult due to the requirement of phase 

separation. Electrical conductivity based methods are 
also generally difficult to apply to multiphase fluid 
mixtures having very lov/ electrical conductivity. The 
capacitance based methods can be applied only to 

20 multiphase fluid mixture where the continuous phase is 
nonconductive. In the case with aqueous slurries, the 
high electrical conductivity of the aqueous phase 
interferes with the dielectric measurement. 

Xn determining the water content of oil/water 

25 emulsion mixtures, prior art systems have a significant 
limitation because of the fact that the electrical 
properties of water -continuous and oil -continuous 
emulsions are quite different even if the water content is 
identical. Prior art systems have also failed to provide 

30 methods or means for determining phase composition in 

fluid mixtures with more than two phases. This is because 
that different sets of phase compositions may result in 
similar conductivity or capacitance measurements. Most of 
prior art systems failed to give accurate measurements 

35 when the concentration of the disperse phase in a two 
phase fluid mixture is low. 

The limitation of conductivity or capacitance 
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based methods is attributed to the limited information 
obtained at a single frequency of excitation alternating 
current (AC) signal. One knovna value of mixture 
conductivity or capacitance is insufficient to determine 
5 the phase composition when both the phase coiqposition and 
the electrical properties of one of the phases in the 
mixture are unknown. 

In certain industrial processes, such as dense 
medium separation of coal and mineral ores and grinding 

10 circuits in mineral processing industry# it is desirable 
to monitor the average particle size of suspended fine 
particles in an aqueous slurry imder the condition of high 
solids concentration. At present there are no simple 
commercially available on-line particle size monitors 

15 capable of this measurement. The conventional method of 
measuring particle size distribution is to remove samples 
from the streams of interest and to perform screen 
analyses on these samples. Hovjever, screen analysis can 
provide a reasonably accurate determination of particle 

20 size distribution above about 45 microns. There are three 
commercially available on-line particle size analysers 
based on ultrasonic attenuation, a scanning laser 
microscope and a reciprocating caliper. However, these 
analysers are not suitable for use in slurry mixtures 

25 where the average particle size is below 45 micron or the 
solids concentration is high or the fluid mediiun is not 
transparent . 

Froth flotation is widely used for concentrating 
minerals, or other valuable constituents, from their ores 

30 or other raw materials. Minerals are separated from 

gangue particles by taking advantage of their differences 
in hydrophobicity . These differences can occur naturally, 
or can be controlled by the addition of a collector 
reagent. Froth flotation generally involves the use of 

35 air injection through a slurry that contains water, 

minerals and gangue particles within a vessel. Dispersed 
air biibbles attract the hydrophobic valuable minerals and 
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carry them upward to the top of the flotation cell, 
whereupon they form a froth bed or froth layer which 
contains and supports pulverised mineral. The froth is 
then scraped or permitted to flow over the lip of the cell 
5 to effect the separation. The thus concentrated mineral 
bearing froth is collected and further processed to 
improve the concentration of desired minerals. The pulp 
may be further processed to recover other valuable 
minerals. 

10 On-line measurement of process parameters is a 

prerequisite for froth flotation process control. 
Whereas some process parameters can be monitored on-line 
with cost effective and reliable measuring devices, the 
effective on-line monitoring and optimal control of froth 

15 flotation processes are still far from being achieved 

because of the strong inertia of the flotation process, a 
still inadequate knowledge of suitable variables for the 
on-line monitoring of the process efficiency and the lack 
of appropriate on-line measurement instrijmentation. 

20 The froth phase in a froth flotation process has 

a nximber of characteristics, including bubble size, 
stability, mobility, solids content and water content. 
The effects of operating conditions such as reagent type^ 
reagent dosage, water chemistry, pulp level, feed flowrate 

25 and aeration rate are reflected in the froth 
characteristics • 

The characteristics of froth layer are related to 
flotation grade and recovery. In view of the difficulty 
in the direct measurement of the froth characteristics, it 

3 0 is desiradDle to use other froth properties that can be 
easily on-line measured as monitoring tools and are 
closely related to flotation grade and recovery. 
SUMMARY OF THE INVENTION 

The present invention provides an alternative 

3 5 method and apparatus for characterising multiphase fluid 
mixtures preferably where the components and mixtures may 
be stationary, moving in bathes or flowing continuously in 
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a conduit/ and more particularly a method and apparatus 
for determining the proportion of each phase constituting 
a multiphase fluid mixture, the type of oil/water emulsion 
mixtures/ the particle size of fine particles in liquid- 
5 solids slurries and the characteristics of bubble froth 
phase* 

In the description hereafter, ^^electrical 
impedance spectrum" refers to the complex plane plot of 
imaginary verses real impedance values for a plurality of 
10 different frecjuencies of energy or in the plotting of 

quantities derived from the real and imaginary impedance 
values • 

According to one aspect of the present invention 
there is provided a method for determining at least one 

15 characteristic of a multiphase fluid including the steps 
of applying alternating energy of a predetermined 
amplitude to a portion of a multiphase fluid and measuring 
the electrical impedance spectrxam across the portion of 
multiphase fluid whereby a characteristic of the 

20 multiphase fluid can be determined from the measured 
electrical impedance spectra w 

It is preferred that the above method is repeated 
for a plurality of different amplitudes of alternating 
energy . 

25 It is preferred that the alternating energy 

includes alternating voltage and alternating current. 

Preferably the electrical impedance spectrum is 
measured across the portion of multiphase fluid for an AC 
voltage of constant amplitude or AC current of constant 

30 amplitude. 

Preferably the alternating energy is applied 
across electrodes in the portion of multiphase fluid. 

The term ^^electrodes'' should be interpreted in 
its broadest sense to include any terminal, wires, or 
35 similar points across which current or voltage can be 
applied to measure the electrical impedance spectriim. 

It is preferred that the electrodes are set a 
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predetermined distance apart and electrical impedance 
measurements are made at the predetermined distance of 
separation between electrodes . 

According to one preferred embodiment of the 
5 invention there is provided an apparatus for 

characterising multiphase fluid laixtures, the apparatus 
including: 

an electrode pair comprising at least one 
conductive path and defining therebetween and thereabout a 
10 sample zone within the multiphase fluid mixture, 

a measuring meajis for measuring characteristics of the 
electrical field formed between the electrode pair, and 

computing means for collecting information from 
the measuring means and converting it to a desired form of 
15 output. 

According to one embodiment of the invention the 
apparatus includes two electrodes. 

Preferably the EIS is measured across the 
electrodes for a constant amplitude of potential 
20 difference (voltage) . 

According to another embodiment the apparatus 
includes three electrodes. 

Preferably the current between adjacent 
electrodes is set at a predetermined amplitude. 
25 It is preferred that the voltage is measured 

across the three electrodes. 

According to a further embodiment of the present 
invention the apparatus includes four electrodes. 

Preferably the four electrodes comprise two pairs 
3 0 of electrodes each adapted to provide a constant magnitude 
of current between the pairs of electrodes. 

Preferably the measuring means is adapted to 
measure the change in voltage between the pairs of 
electrodes. 

35 It is preferred that the apparatus includes a 

configuration of electrodes in which either a current of 
constant amplitude is applied across the electrodes, a 
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voltage of constant amplitude is applied across the 
electrodes or a combination of current of constant 
amplitude and voltage of constant amplitude is applied 
across respective adjacent electrodes. 
5 It is preferred that a current of cpnstant 

amplitude is provided by grouping a pair of electrodes. 

According to another embodiment it is preferred 
that electrical impedance spectra are measured for 
different values of constant voltage or alternatively 
10 different values of constant current. 

According to one aspect of the present invention 
there is provided a method for characterising a multiphase 
fluid mixture, the method including the steps of 

applying an alternating current or voltage to the 
15 electrodes located in the multiphase fluid mixture, 

measuring the electrical impedance spectrum 
across the electrodes at one or a few selected amplitudes 
of excitation signal; 

transforming the measured electrical impedance 

2 0 spectrum or spectra into a few indicator quantities using 

feature extraction algorithms; and 

determining, from the indicator quantities of the 
measured electrical impedance spectrum, at least one 
characteristic of at least one phase constituent in the 
25 mixture using a predictive mathematical model. 

According to one embodiment, the invention 
involves measuring the real and imaginary parts of the 
impedance over a frequency range of O.lHz to IMHz. Real 
and imaginary impedance values preferably include real and 

3 0 imaginary components of mathematically related parameters 

such as impedance, admittance, modulus and dielectric 
permittivity, etc. The impedance sensing means may be 
configured in two or three or four electrodes. 

It is preferred that the method includes the step 
35 of measuring temperature and pH value of the multiphase 
fluid mixture at each measurement frequency or measured 
impedance spectra. 
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The multiphase fluid mixtures preferably include 
matter such as gas, solid, liquid or different 
combinations of the above. 

It is preferred that the method includes the step 
5 of transforming the measured electrical impedance spectnun 
into a few indicator quantities- The step of ' transforming 
further includes obtaining data in the form of average 
number of good readings, calculating a smoothed electrical 
impedance spectrum using 2D data smoothing algorithms, 

10 such as locally weighted regression, and scaling indicator 
quantities and temperature and pH readings into suitable 
value ranges. The indicator quantities of the electrical 
impedance spectrum may further include the real and 
imaginazry impedance values at a number of selected 

15 frequencies, first and second derivatives of the spectrtim 
at a number of selected frequencies, average values of 
imaginary impedance component over a selective range of 
real impedance, the parameters of a mathematical model for 
the representation of the impedance spectrum, and the 

20 latent variables, summarising the information contained in 
the original impedance spectrum, calculated from 
multivariate statistical methods, such as principal 
component analysis (PCA) and partial least-squares (PLS) . 
The method preferably includes calculating 

25 indicator quantities by fitting the electrical impedance 
spectrum to a mathematical model of the impedance 
spectrum. The mathematical model may further include an 
electrical equivalent circuit model and empirical 
regression equations. 

3 0 Preferably the method includes the step of 

determining, from the indicator quantities combined with 
temperature and pH, at least one characteristic of at 
least one phase constituent in the mixttire using a 
predictive mathematical model. The characteristics 

3 5 further . include the proportion of each phase constituting 
a multiphase fluid mixture, the type of oil/water emulsion 
mixtures, the particle size of fine particles in liquid- 
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solids slurries and the characteristics of bubble froth 
phase. The predictive mathematical model further include a 
trained artificial neural network and a multivariate 
regression model. 
5 It is preferred that the method includes training 

and validating an artificial neural network with a number 
of indicator quantities with known characteristics of 
multiphase fluid mixtures. The method further includes 
calculating parameters in a predictive mathematical model 
10 using a number of indicator quantities with known 
characteristics of multiphase fluid mixtures. 

It is preferred that the method includes the step 
of analysing the impedance spectrum using pattern matching 
algorithms to determine whether characteristics of a 
15 bubble froth phase loaded with particles are favorable or 
not in terms of the grade and yield of the flotation 
concentrate . 

The method preferably includes the step of 
analysing the impedance spectrum using pattern matching 
2 0 algorithms to determine whether an oil/water emulsion is 
oil continuous or water continuous type. 

According to a further aspect of the present 
invention there is provided a method of analysing 
extraneous matter in a fluid including the steps of 
25 receiving impedance data, being data including real and 
imjaginary impedance values measured across electrodes 
located in a fluid, recording the impedance spectrum at a 
plurality of time intervals, calculating indicator 
quantities of the impedance spectrxam for the received 
30 impedance spectrum data, comparing indicator quantities of 
the impedance spectrum with reference indicator quantities 
and determining at least one characteristic of at least 
one phase constituent in the mixture from the comparing 
steps . 

35 The method may include the step of determining 

characteristics of multiphase fluid mixtures in the forms 
of numerical value or qualitative index. 
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It should be noted that reference to electrical 
impedance spectruia refers to EIS (Electrical Impedance 
Spectrum) . 

It is to be understood that, if any prior art 
5 p\iblication is re'ferred to herein, such reference does not 
constitute an admission that the publication ^orms a part 
of the common general knowledge in the art, in Australia 
or in amy other coxantry. 

The words comprising., having, including" should 
10 be interpreted in an inclusive sense, meaning that 
additional features may also be added. 
BRIEF DESCRIPTION OF DRAWINGS 

Preferred embodiments of the present invention 
will now be described by way of example only with 
15 reference to the accompanying drawings in which: 

Figure 1 is a block diagram of the measurement 
apparatus for characterising multiphase fluid mixtures 
according to a first embodiment of the present invention; 

Figures 2A, 2B, 2C and 2D are schematic diagrams 
20 illustrating possible designs of electrode pairs useful 
for the measurement of electrical impedance spectrum; 

Figures 3A and 3B show graphical 
representations of electrical impedance spectra of 
liquid- solids slurries v/ith different compositions; 
25 Figures 4A and 4B show graphical 

representations of electrical impedance spectra of sugar 
syrup with different sugar crystal contents; 

Figures 5A and 5B show graphical 
representations of electrical impedance spectra of water- 
30 oil emulsions; 

Figure 6 shows graphical representations of 
electrical impedance spectra of liquid- solids slurries 
with different particle sizes; 

Figures 7 shows graphical representations of 
35 electrical impedance spectra of froth (or foaming) phase 
with different characteristics; 

Figure 8 shows a schematic diagram of a method 
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for characterising multiphase fluid mixtures according to 
one embodiment of the present invention; and 

Figure 9 shows a PCA neural network to project 
the data from D to M dimensions. 
5 ' It should be \mder stood that the em]bodiments of 

the invention described hereinafter with reference to the 
drawings refer to specific electrode configurations where 
the electrode type, nxamber of electrodes and distance 
between electrodes remains fixed. The invention also 
10 covers other embodiments where different n-umbers and 
distances between electrodes are provided as well as 
different types of electrodes.. In these other embodiments 
values of electrical impedance would be different to those 
exemplified in the preferred embodiments. 
15 DETAILED DESCRIPTION 

According to one embodiment of the present 
invention the electrical impedance spectrum of a 
multiphase fluid mixture was measured over a wide range of 
frequencies to identify characteristic parameters of 
2 0 interest in a multiphase fluid mixture. In addition the 
inventors noted the dependence of electrical properties of 
constituents in multiphase fluid mixtures upon excitation 
by an AC signal varies. It was therefore considered that 
the electrical impedeoice spectrum of a multiphase fluid 
25 mixture measured over a wide range of frequency may 

contain sufficient information for deducing characteristic 
parameters of interest. 

The inventors also realised that electrical and 
dielectric properties of solids-liquid suspensions depend 
30 on not only the phase composition but also the particle 

size of solids. When an AC current is passing through a 
suspension, the surface charge and the associated 
electrical double layer of particles tend to cause a phase 
shift of the AC current in a certain range of frequency 
35 due to charge relaxation processes on the surface. For a 
given volume fraction of suspended particles, the smaller 
the particle size, the higher the amount of the surface 
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charge- Since the phase shift is proportional to the 
amount of the surface charge, small particles will cause 
a higher phase shift than large particles. It is, 
therefore, possible. to calculate the particle size from 
5 measured real and imaginary parts of electrical impedance 
over a wide range of frequency. 

Furthermore the inventors discovered that the 
electrical and dielectric properties of components in a 
froth phase are different from each other. 

10 From the viewpoint of electrical behaviour of the 

froth phase, the inter-b\abble lamellae containing water 
and solids can be regarded as a complex network of 
electrical conductance, inductance and capacitance. The 
structure of this network would be sensitive to changes in 

15 the froth structure and characteristics. Therefore, the 

measurement of the electrical impedance of the froth phase 
over a wide range of frequencies would probe into the 
froth structure and/or characteristics. 

As shown in Figure 1 an apparatus for 

20 characterising multiphase fluid mixtures consists of a 

pair of fluid measurement electrodes 11a and lib immersed 
in multiphase fluid mixtures 12 ^ a temperature sensor 13, 
cm EIS and temperature measurement unit 14, a computing 
unit 15 and an output unit 16. 

25 Referring to Figures 2A-2D/ the measurement 

electrodes can be mounted on the inner surface of a 
conduit or vessel wall in the forms of tapped rods 20a and 
20b or rings 23a and 23b, Alternatively the electrodes 
can be mo\inted on a non- conductive rod 24 in the form of 

30 dots 25a and 25b, or on a non -conductive spacer 26 in the 
form of plates 27a and 27b with any suitable shapes. 
Instead of the plate type of electrodes 27a and 27b, one 
electrode may be a rod electrode surro\inded coaucially by 
smother cylindrical electrode. In the cases with a 

35 conductive conduit or vessel wall, the electrodes 2 0a and 
20b have to be insulated with the conductive wall 21 and 
a non- conductive layer 22 has to be applied to cover the 
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inner surface of the conductive wall 21. The material for 
non- conductive layer 22 includes certain ceramics, casting 
basil, plastics and other suitable materials. 

An electrical impedance spectrum and temperature 
5 measurement unit 14 is connected to each of electrodes as 
well as to the computing unit 15. The measurkcient unit 14 
sends and receives signals to or from the computing unit 
15 through electrical, optical, electromagnetic wireless 
or other type signals. The output unit 16 preferably is a 

10 visual displayer, e.g. LCD, for displaying the results 
provided by the computing \anit 15. 

The measurement unit 14 preferably includes a 
signal generation module for generating AC signals at 
specified amplitude and frequencies, a measuring module 

15 for measuring the amplitude and phase angle of AC signals, 
a temperature measurement circuit, self calibration and 
diagnosis circuits and an embedded microprocessor for 
controlling signal generation and measuring module and 
sending and receiving signals to or from the computing 

20 vmit 15. 

The computing unit 15 preferably includes means 
for outputting control variables or commands to the 
measurement unit 14, means for receiving and recording 
measured temperature, real and imaginary impedance values 

25 for a plurality of different frequencies, means for 
checking the validity of received data, means for 
scaling the received data into a suitable value range, 
means for calculating indicator quantities from the 
measured EIS, means for clustering the data into data 

30 patterns and means for determining at least one of 

characteristics of multiphase fluid mixtures from the 
indicator quantities . 

By measuring the electrical impedance spectrixm 
across the fluid measurement electrodes 11a and lib, 

35 information about characteristics of multiphase fluid 

mixtures can be identified. The multiphase fluid mixtures 
preferably include matter such as gas, solid, liquid or 
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different combinations of the above. 

For example as sbovm in Figures 3A and 3B the 
effects of phase composition on electrical impedance 
spectruun of slurries can be ascertained. Electrical 
5 ' impedance spectra for water only, water slurry containing 
20% (by volume) sands, water slurry containing 12% 
magnetite and water slurry containing 10% magnetite and 
20% sands are represented by 31, 32, 33, 34, respectively. 
Figure 3A shows the spectra for the frequency range of 0.5 

10 Hz to IMHz^ In order to emphasise the effects of phase 

composition on EIS, the same spectra only in the frequency 
range of 100 Hz to IMHz is shown in Figure 3B. Figures 
3A cuxd 3B clearly indicate that the EXS is sensitive to 
the changes of phase composition of aqueous slurry 

15 mixtures. It is this sensitivity that provides the basis 
for the present invention. It can be also seen that the 
effect of the presence of magnetite on the spectrum is 
substantially different from that of sand. The presence 
of magnetite can cause a peak in the high frequency range 

20 of the spectrum, but sand cannot. The ability of the 
apparatus in the present invention to distinguish the 
relative composition of different dispersed phases is 
based on their different effects on the spectnun* The 
measurement of the spectra as shown in Figures 3 A and 3B 

25 can be repeated at a few different aiirplitudes of the 

excitation signal and the determination of the amplitude 
dependence of the spectra would allow the further 
differentiation of factors causing the changes of the EIS. 

In Figures 4 A and 4B EIS is produced for sugar 

30. syrup having different crystal contents. In the case 
with white sugar (see Figure 4A) , the EIS for the 
unsaturated symip 41 (containing 20% by volume water and 
80% saturated syrup) is significantly different from that 
for saturated syrup 42. By adding 10% by weight white 

35 sugar crystals into the saturated syrup, a spectrum 43 is 
produced. Electrical impedance spectra for unsaturated 
raw sugar syrup, saturated raw sugar syrup and saturated 
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raw sugar syrup with 10% (by weight) raw sugar crystals 
are represented by 44, 45 and 46, respectively. By 
comparing the EIS for raw sugar in Figure 4B with those 
for white sugar in Figure 4A, it can be seen that the 
5 spectra for raw sugar syrup have a lower value of real 
impedance and different spectrxim patterns from those for 
white sugar. This is due to the higher concentration of 
soluble impurity in raw sugar. As shown in Figures 4A and 
43, EIS can detect not only the crystal content of mother 

10 liquor but also the purity of mother liquor. 

Figures 5A and 5B show examples of EIS curves for 
water-oil emulsions • Curve 51 shows the EIS for the 
water-in-oil emulsion with 25% (by volume) water whereas 
curve 52 shows the EIS for the emulsion with 50% water. 

15 Figure 5B shows EIS for oil-in-water emulsions, and the 
spectra for 50% and 7 5% water are represented by 53 and 
54, respectively. It can be seen that the EIS pattern for 
water- in-oil emulsions is different from that for oil-in- 
water emulsions. This difference will provide a basis for 

2 0 identifying emulsion type using EIS. 

Figure 6 shows the EIS change of slurry with 
particle size under the same volumetric concentration of 
solids. It can be seen that for this particular 
particles the EIS for 30 |im particle sise 61 is 

25 significantly changed to curve 62 when the particle size 
is reduced to 20 pm. It should be pointed out that EIS 
is not sensitive^ to particle size change when the size is 
higher than 50 pm. However, in situations where the 
particle size is smaller than 50 |.im, it is possible to 

30 monitor the particle size by observing the change of EIS 
curve i 

The froth phase in bubble flotation processes is 
a special type of multiphase fluid mixtures, in which the 
electrical and dielectric properties of components are 
35 different from each other. For example, the conductivity 
of water is several orders of magnitude higher than that 
for mineral particles. From the viewpoint of electrical 
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behaviour of the froth phase , the inter-bubble lamellae 
containing water and solids can be regarded as a complex 
network of electrical conductance, inductance and 
capacitance. The structure of this network is sensitive 
5 to changes in the operating conditions of bxibble' flotation 
processes, and hence the effects of operating' conditions, 
such as reagent dosages, feed flowrate and froth depth, on 
the flotation performance are reflected on the measured 
EIS. Therefore, the measurement of the electrical 

10 impedance of the froth phase over a wide range of 
excitation signal frequency would probe into the 
performance of flotation processes • 

Figure 7 shows electrical impedance spectra 
measured in the froth phase of bubble flotation processes 

15 of one fine coal under various operating conditions. The 
spectra for 7 8%, 74% and 68% flotation yield are 
represented by 71, 72 and 73, respectively. It can be 
seen from the figure that the EIS spectra is closely 
correlated with the product yield. For this particular 

20 coal, the spectrum 71 is favourable in term of product 
yield. This favourable spectrum pattern can be 
conveniently used as the objective function for optimising 
operating conditions. In the bubble flotation of other 
materials, such as minerals, the pattern of EIS of the 

25 froth phase may be different from that shown in Figure 7. 
However, the favourable pattern of EIS and associated 
operating conditions still can be identified using EIS as 
long as the flotation performance is sensitive to the 
changes in operating conditions • 

30 Examples presented in Figures 3 to 7 clearly 

demonstrate that the electrical impedance spectr\im can 
provide sufficient information regarding to the 
characteristics of multiphase fluid mixtures. In order to 
use these information for the on-line estimation of the 

35 characteristics of multiphase fluid mixtures, a 

mathematical or other type of relationship between the 
EIS and its corresponding characteristics of multiphase 


BNSDOCID: <W O 20Q4077036A1 1 > 


wo 20()4/077()3r> 


- 17 - 


PCT/AU2004/00()187 


fluid mixtures is required. Among the various approaches 
for describing and modeling phenomena that are too complex 
for analytical methods or empirical rules, artificial 
intelligent data analysis techniques, particularly the 
5 artificial neural network (ANN) ' have shovm great potential 
as an effective method for identifying or mapping complex 
non- linear relations without rec[uiring specific laiowledge 
of the model structure. Artificial neural network 
techniques are very efficient in computation due to the 
10 feedforward nature and also have higher tolerance to 
errors in the input data set than other parameter 
estimation approaches. Hence, a multiplayer perceptron 
artificial neural network (MLP-ANN) is a preferred but not 
an exclusive approach in the present invention to estimate 
15 characteristics of interest from the measured EIS of 

multiphase fluid mixtures. Other approaches, such as 
multivariate regression and ANN based on fuzzy logic are 
also useful in correlating the measured EIS with the 
characteristics of multiphase fluid mixtures. 
20 Based on observations derived from EIS 

measurements taken using the aforementioned apparatus it 
is possible to employ an automated procedure to identify 
characteristics of multiphase fluid mixtures. Figure 8 
is a flow chart showing a method for implementing this 
25 automated procedure. 

As illustrated in Figure 8, when the power is ON, 
the measurement unit 14 makes a diagnosis of itself and 
becomes initialised in step 80. Computing unit 15 then 
sends control variables to the measurement \init 14 in step 
30 81. Control variables include the amplitude of AC signal 
generated by 14, frequency range, number of measurement 
points in the frequency range, and the like. 

Once the measurement unit 14 receives control 
variables the electrical impedance spectrum, temperature 
3 5 and optionally pH are measured and recorded in step 82 

using the aforementioned apparatus. It is preferred that 
the measurement of EIS, temperature and pH in step 82 are 
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repeated several times in a short period of time and their 
average values are used for further processing. If the 
data is valid for a particular application as referenced 
by step 83 the computing unit is able to activate a data 
5 processor so as to scale the data. into a suitable range of* 
values as referenced by step 84. Alternatively if the 
data is not valid an alarm signal is provided to a display 
to notify an observer that the invalid data occurs and 
the measurement and recording step 82 is repeated. 

10 After the data has been scaled into a suitable 

range the computing unit 15 is programmed to calculate the 
indicator quantities from the scaled EIS data as 
referenced by step 85. Then a software program performs 
a classification analysis of data pattern in step 86 to 

15 identify whether the EIS data pattern is unseen in the 

training stage of an artificial neutal network (ANN) or in 
the development stage of a multivariate regression model. 
If the answer is yes an alarm signal is provided to a 
display to notify an observer that the new data pattern 

20 occurs amd the measurement and recording step 82 is 

repeated. If the new data pattern repeatedly occurs ^ the 
computing unit 15 is programmed to retrain an ANN model or 
refit a multivariate regression model using a data set 
including the new data pattern. ^Alternatively if the 

25 data pattern is not a new one an output of at least one of 
the characteristics of the multiphase fluid mixtures is 
produced by the computing unit in the step 87. If there 
is no manual interruption then the measurement and 
recording step of item 82 is repeated. 

3 0 In the data validation step, as referenced by 

item 83, data with a low precision, values close to pre- 
specified limits and significant noise are discarded to 
control the data quality for further processing. In ANN 
and multivariate analysis it is mandatory to scale the 

35 measured EIS and other data before the main business of 
analysis begins. This is because the measured EIS and 
other variables have different units and magnitude of 
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10 


values. Scaling methods useful in the present invention 
include column centring, standardisation and range 
scaling. Range scaling cause the values to fall between 0 
to 1 or -1 to 1. These scaling methods are applied only 
to columns (i.e. data points at a same frequ.ency from 
different measurements) . 

In order to capture all important frequency and 
signal amplitude dependent information, a number of 
frequency points are usually used in the measurement of 
EIS and the measurement is repeated with a few different 
amplitude of excitation signal. In applying mathematical 
approaches, such as artificial neural network and 
multivariate regression, to predict characteristics from 
measured EIS, the use of all data points in a spectrum 
15 will result in a very large dimension of input. An 

unnecessary large dimension of input variables will have 
adverse effects. For a fixed n\amber of training data 
patterns, with the increase of input variables it becomes 
more sparse in the multi -dimensional space, and therefore 
20 degrades the learning performance. The generality of the 
trained ANN model may also be reduced due to inclusion of 
irrelevant or unimportant input variables. Apart from 
irrelevant and unin^ortant variables that cause large 
dimension of input variables, there may be correlation's 
25 between EIS data points measured at frequencies close to 
each other. Correlated inputs make the model more 
sensitive to the statistical peculiarities of the 
particular data sample, and they accentuate the 
overfitting problem and limit generalisation. Therefore, 
30 it is an important step in the present invention to 

calculate or extract from EIS indicator quantities with a 
much less n\imber of variables but retaining sufficient 
information of the original spectrum. 

The indicator quantities of the electrical 
35 impedance spectrum may include the real and imaginary 
impedance values at a ntamber of selected frequencies, 
first and second derivatives of the spectrum at a number 
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of selected frequencies, average values of imaginary 
impedance component over a selective range of real 
impedance, the parameters of a mathematical model for the 
representation of the impedance spectrum, and the latent 
5 variables or principal components, sximmarising the 

information contained in the original impedance spectriim, 
calculated from multivariate statistical methods, such as 
principal component analysis (PCA) and partial least- 
squares (PLS) • 

10 The method preferably includes calculating 

indicator quantities by fitting the electrical impedance 
spectrum to a mathematical model of the impedance » 
spectrum. The mathematical model may further includes an 
electrical equivalent circuit model and empirical 

15 regression equations. 

It is preferred that PCA Implemented using an 
artificial neural network, as shown in Figure 9, adapted 
with Hebbian learning or similar rules is used for 
calculating indicator quantities for the robustness. 

20 There are well-known algorithms that analytically computer 
PCA, but they have to solve matrix equations associated 
with singular value decomposition. When the matrices are 
ill-conditioned, the numerical solutions fail, while PCA 
neural networks provide more robust solutions. 

25 Xf an indicator quantity data set has a pattern 

which has not been seen in tlie training stage of an ANN 
model or In the development stage of a multivariate 
regression model, the output of these model taking the 
data set as input will be erroneous. Therefore, it is 

3 0 necessary to check whether the pattern of a new data set 
Is new. The classification of data patterns can be- 
performed using ANN based approaches, such as unsupervised 
Bayeslan clustering system, or the data reconstruction 
approach associated with PCA. If the indicator quantities 

3 5 calculated from a measured EIS (x) is represented by 
vector y, the reconstructed EIS is given by 
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where W is the weight matrix. 

If the difference between the reconstructed EIS, 
X' and the measured EIS, x is larger than a specified 
threshold, the data pattern of the measured EIS can be 
considered as a new one. If this new pattern, repeatedly 
occurs, it will become necessary to retrain an ANN model 
or refit a multivariate regression model. 

In the prediction step, as referenced by item 87, 
the ANN model can be replaced by a multivariate regression 
model, a pattern matching algorithm or even a lookup 
table. When a lookup table is used in step 87, the 
con^arison between indicator (juantities of a measured EIS 
with reference indicator quantities will be used to 
determine characteristics of multiphase fluid mixtures. 
3^5 The output from an trained ANN model or other 

types of relationship, such as lookup tables can be 
numerical values or qualitative indices, such as 
classification index. 


10 
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CliAIMS 

1. A method for determining at least one 
characteristic of a multiphase fluid including the steps 
of applying alternating energy of a predetermined 
5 amplitude to a portion of a multiphase fluid ..and measuring 
the electrical impedance spectarum across the portion of 
multiphase fluid whereby a characteristic of the 
multiphase fluid can be determined from the measured 
electri.cal impedance spectra. 

10 2. The method as claimed in claim 1 including 

the step of applying alternating energy of a different 
predetermined amplitude to the portion of the multiphase 
fluid and measuring the electrical impedance spectrum 
across the portion of the multiphase fluid whereby the 

15 characteristic of the multiphase fluid can be determined 
from the measured electrical impedance spectra. 

3 . The method as claimed in claim 1 Including 
the steps of applying a plurality of alternating currents 
of different constant amplitude to the portion of the 

20 multiphase fluid and measuring the electrical impedance 
spectrum across the portion of multiphase fluid whereby 
the characteristic of the multiphase fluid can be 
determined from the measured electrical impedance spectra 
for each qf the constant amplitudes. 

25 4. The method as claimed in claim 1 including 

the steps of applying a plurality of alternating voltages 
of different constant amplitude to the portion of 
multiphase fluid and measuring the electrical impedance 
spectrum across the portion of multiphase fluid whereby 

30 the characteristic of the multiphase fluid can be 

determined from the measured electrical impedance spectra 
for each amplitude of constant voltage. 

5. The method as claimed in claim 1 including 
the step of providing at least two electrodes and applying 

35 alternating energy of the predetermined amplitude to the 
portion of multiphase fluid between the electrodes. 

6. The method as claimed in claim 1 including 
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the step of providing three or more electrodes to apply 
the alternating energy across. 

7 . The method as claimed in claim 1 including 
the step of transforming the measured electrical impedance 
5 spectrum into a plurality of indicator values, using a 
feature extraction algorithm. 

8 • The method as claimed in claim 7 including 
the step of determining from the indicator values at least 
one characteristic at least one phase constituent in the 

10 multiphase fluid. 

9 . The method as claimed in claim 8 including 
the step of using predictive mathematical modelling to 
determine one or more characteristics of the multiphase 
constituents . 

15 10 . The method as claimed in claim 1 including 

the step of measuring the real and imaginary parts of the 
impedance over 0 . IHz to 1 MHz . 

11. The method as claimed in claim 10 including 
the step of measuring temperature and pH values of the 

20 multiphase fluid at each measured spectrum. 

12 . The method as claimed in claim 11 wherein 
the transforming step includes the step of calculating a 
smoothed electrical impedsmce spectrum using a smoothing 
algorithm. 

25 13. The method as claimed in claim 12 wherein 

the step of transforming includes calculating indicator 
values by fitting the electrical impedance spectrum to a 
matheiaatical model of the impedance spectrum. 

14 . The method as claimed in claim 13 wherein 
30 the mathematical model includes an electrical equivalent 

circuit model and emperical regression equations. 

15. The method as claimed in claim 14 wherein 
the step of determining includes using a predictive 
mathematical model . 

35 16. The method as claimed in claim 15 wherein 

the predictive mathematical model includes a trained 
artificial neural network and a multi variate regression 
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model • 

17 . The method as claimed in claim 16 including 
the step of training and validating an artificial neural 
network with a number of indicator quantities with known 

5 characteristics of multiphase fluid mixtures.., 

18. The method as claimed in claim' 17 including 
the step of analysing the impedance spectrum using pattern 
matching algorithms to determine whether characteristics 
of a bubble froth phase loaded with particles are 

10 favourable or not in terms of the grade and yield of the 
flotation concentrate. 

19. The method as claimed in claim 17 including 
the step of analysing the impedance spectrum using the 
pattern matching algorithms to determine whether an 

15 oil/water emulsion is of an oil continuous or water 
continuous type . 

20. An apparatus for characterising multiphase 
fluid mixtures, the apparatus comprising at least one pair 
of electrodes for meausuring at least one characteristic 

20 of a sample zone of a multiphase fluid located 

therebetween, a field generation means for generating an 
electrical field between the electrodes, a measuring means 
for measuring at least one characteristic of the 
electrical field formed between the electrodes and a data 

25 processing means for collecting data from the measuring 
means, processing the data and outputting real axid 
imaginary impedance data for a constant amplitude of 
voltage or current generated by the field generation 
means . 

30 21. The apparatus as claimed in claim 20 

including three electrodes configured to have a constant 
voltage between first and second electrodes and a constant 
voltage between second and third electrodes. 

22. The apparatus as claimed in claim 20 

35 including four electrodes configured to have a constant 
current applied by first and second electrodes and by 
third and fourth electrodes. 
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23. An apparatus as claimed in claim 20 
including a plurality of pairs of electrodes configured to 
have a constant current produced by some electrodes and a 
constant voltage across other electrodes. 
5' 24. The apparatus as claimed in cl,aim 20 

wherein the field generation means is configured to 
generate frequencies between O.lHz and 1 MHz. 

25- The apparatus as claimed in claim 24 
including a temperature sensor for measuring the 
10 temperature of multiphase fluid within the sample zone and 
a pH value sensor for sensing the pH value of multiphase 
fluid within the sample zone- 

26. A method of analysing extraneous matter in 
a fluid including the steps of receiving impedance data, 

15 being data including real and imaginary impedance values 

measured across electrodes located in the fluid, recording 
the iinpedance spectrum at a plurality of time intervals at 
a predetermined amplitude of energy applied across the 
electrodes, calculating indicator quantities of the 

20 impedance spectrum for the received impedance spectrum 
data, comparing indicator quantities of the impedance 
spectrum with reference indicator quantities and 
determining at least one characteristic of at least one 
phase constituent in the multiphase fluid for the 

25 comparing steps. 

27. The method as claimed in claim 26 wherein 
the indicator quantities include the minimxjm nxamber of 
quantities required to model the original electrical 
impedance spectrum. 

30 28. The method as claimed in claim 27 wherein 

the indicator quantities include one or more of the real 
and imaginary impedance values at a number of selected 
frequencies, first and second derivatives of the spectrtam 
at a number of selected frequencies, average values of 

35 imaginary impedance component over a selective range of 

real impedance, the parameters of a mathematical model for 
the representation of the impedance spectrum and the 
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latent variables or principle components summarising the 
information contained in the original impedance spectrum 
calculated from at least one multivariate statistical 
method. 

5 29. The method as claimed in claim^ 28 wherein 

the multivariate statistical method includes a principle 
component analysis using an artificial neural network. 

30. The method as claimed in claim 29 including 
a checking means to check whether a pattern of a new data 

10 set of indicator quantity data fits a predetermined 
pattern. 

31. The method as claimed in claim 30 wherein 
the checking means includes the step of representing 
indicator quantities calculated from a measured electrical 

15 impedance spectrum EIS (x) by vector y and calculating the 
reconstructed EIS given by x' = W^y where W is the weight 
matrix, whereby if for the difference between the 
reconstructed EIS, x' and the measured EIS, x is larger 
than a specified threshold, the data pattern of the 

20 measured EIS is recorded as a new one. 
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